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1. INTRODUCTION 

Renewable energy refers to the energy derived from sources that are naturally replenished by nature 
like the sunlight for solar energy, the wind that powers turbines, the water for hydroelectric energy, and 
geothermal heat. These forms of energy generation are widely accepted as being vital for the continuance of 
life on planet earth [1, 2]. As the world`s population grows larger, the energy requirements also increase. There 
was a 50% increase in human energy consumption from 1980-2003, which was caused by a growth of 42% in 
the global population, especially in developing countries. From this fact, it can be estimated that the demand 
for energy will increase by 57% by 2025. Indonesia is a country with the potential to generate large amounts 
of geothermal energy because it is located in the ring of fire and there are numerous active volcanoes on 
the tectonic plate that runs from Sumatra to Java, and from Bali and Maluku up to Sangihe islands. 
Being located near volcanoes provides high enthalpy values in terms of the geothermal energy potential. 
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The geothermal resource refers to the amount of geothermal heat that can be generated if specific technological 
and economic conditions are met [3]. The Indonesian government stated at the 2000 World Geothermal 
congress that there were 276 locations that could be used for geothermal energy production. The locations were 
in Sumatra, Java, Bali, Kalimantan, Nusa Tenggara, Maluku, Sulawesi, Halmahera and Papua, and they 
estimated that the total potential energy that could be produced was 27.354 MWe (Mega Watt equivalent). 
More recent estimates of the geothermal potential are closer to 30.000 MWe. However, at present only 196 
MW is generated, which represents around 40% of the total potential. 

Satellite remote sensing (RS) offers a wide variety of image data with different characteristics like 
temporal, spatial, radiometric and spectral resolution. These features are useful because they match 
the requirements for monitoring environmental conditions and natural resources [3, 4]. Nowadays, large 
amounts of data are provided by satellite sensors like Landsat, IKONOS, SPOT-5, WordView-2, QuickBird, 
GeoEye, and Orbview. In geological application, RS satellite images assist in conducting cost-effective 
geothermal exploration and in identifying the geothermal potential by ascertaining the land surface temperature 
(LST) through thermal infrared remote sensing (TIR). Now that more emphasis and resources are 
directed towards the search for renewable energy sources, TIR is regarded as vital when exploring geothermal 
resources [5]. One of the satellite sensors that TIR images can be obtained from is Landsat’s Thermal Infrared 
Sensor (TIRS), which contains a wealth of information about the surface characteristics of the earth over 
the past 40 years. In particular, the thermal emission in band number ten on the sensor (10.60—11.19 um) shows 
the changes in temperature observation capabilities. The images help to determine the thermal characteristics 
of an area and to detect temperature anomalies. The anomalies that indicate a high potential of geothermal 
resources are the areas with a high LST. TIR is acknowledged as an effective method for obtaining quantitative 
and qualitative LST data on regional and global scales [5, 6], and LST is regarded as a key factor in the study 
of geological surface activity [7]. Landsat TIR data can be used for monitoring environmental conditions, 
spatial decision making and monitoring water consumption. Also, it can be used for investigating urban 
microclimates, mapping sensible heat flux, monitoring volcanoes, and observing fire depleted flora by using 
burnt area mapping. Observing LST provides vital data on the climate and physical characteristics of 
the planets surface. This information is used in studies on global warming, evapotranspiration, hydrological 
cycles, vegetation monitoring, urban climatic conditions, and the urban heat island effect [8]. LST is becoming 
more important in environmental studies and different methodologies are being implemented to observe 
LST from space [9]. LST is also mentioned in a number of studies on ecology, hydrology, and global 
change [10, 11]. 

LST specifically refers to the temperature of the Earth's surface, and it is obtained from either satellite 
imagery or by taking direct measurements. Satellite imagery helps to monitor the variation in the earth's 
resources and to observe the physical processes of surface energy and water balance at a local and global scale. 
There are several algorithms for estimating LST from satellite TIR measurements which are based on varying 
assumptions [12, 13]. The resulting LST measurements are crucial for evaporation modelling, climate 
modelling, and radiative transfer modelling [14]. 

This research analyses the land surface temperature to identify potential geothermal resources based 
on remote sensing using Landsat 8 TIRS imagery from satellite data. Since it was provided thermal infrared 
bands (Band 10 and 11) for estimating brightness surface. Single channel/Mono-window algorithms were used 
to generate the LST maps. And data set was combined with a digital elevation model (DEM) to identify 
the potential geothermal energy based on the variation in surface temperature. TIR remote sensing is 
an excellent method for obtaining LST data because the sensor detects most of the energy that emits from 
the surface of the land [15, 16]. 


2. RESEARCH METHOD 

The research area in this case study is located West Sumatra in Indonesia with the geographic 
coordinates of 1°00’'S 100°30’E. The majority of West Sumatra consists of the highlands of the Bukit Barisan 
mountain range which runs through the centre of the province from north to south. There are lowland areas 
surrounding the mountains and the west coast of the province faces the Indian Ocean. West Sumatra borders 
the provinces of North Sumatra and Riau in the north, as well as Jambi and Bengkulu to the south. This research 
focuses on the potential in 5 districts of West Sumatra, they are Pasaman, Agam, Solok, Tanah Datar, and Lima 
Puluh Kota as shon in Figure 1. 

Landsat 8 carries two instruments: the first is TIRS (thermal infrared sensor) which includes band 10 
and 11 and OLI (operational land imaginer) sensor Band (1-9), as shown in Table 1. The instruments collect 
image data from nine shortwave spectral bands (OLI1~ OLI9) over a 185 km swath with a 30 m spatial 
resolution for all bands except for the 15 m panchromatic band (OLI8); TIRS provides image data for both 
thermal bands (TIRS10, TIRS11) with a resolution of 100 m. The data for this study was taken on 01 July 2015 
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from Landsat 8 in a scenel (path 127, row 60) and scene 2 (path 127, row 61). Panchromatic band 8 had 
a spatial resolution of 15 m. Universal Transverse Mercator (UTM) projection with zone 43 North and WGS84 
datum were used for the datasets. Landsat TIR data can be used for observing water consumption, 
environmental conditions, and spatial decision-making. This type of data is also useful for examining 
microclimates in urban areas, mapping sensible heat flux, observing volcanoes, and monitoring fire damaged 
flora through burnt area mapping [17]. Tables and figures are presented center, as shown below and cited in 
the manuscript. Mono window algorithm used in this research to examine land surface temperature with 
Landsat-8 OLI imagery involved seven stages: pre-processing (Landsat-8 OLI imagery, cloud removal, 
ortho-rectification and mosaicking, colour balancing), and processing (processing temperature, extraction land 
surface temperature, top of atmosphere (ToA) spectral radiance, brightness temperature, estimation of 
emissivity, calculation NDVI and land surface temperature estimation), see in Figure 2. 
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Figure 1. Location of study area map in West Sumatra 


Table 1. The Spectral range and spatial resolution of the Landsat 8 OLI and TIRS bands [18] 


Sensor Band Spectral Range (um) Spatial Resolution (m) 
OLI 1 0.43-0.45 30 

2 0.45-0.51 30 
3 0.53-0.59 30 
4 0.64-0.67 30 
5 0.85-0.88 30 
6 1.57-1.65 30 
7 2.11-2.29 20 
8 0.50-0.68 15 
9 1.36-1.38 30 

TIRS 10 10.60-11.19 100 
11 11.50-12.51 100 
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2.1. Pre-processing of Landsat-8 OLI imagery 

The Landsat-8 OLI images were in tagged image file format (TIFF) and the data type was integer. 
During image pre-processing, cloud removal is a common practice because areas that are of interest are 
sometimes obscured by clouds, which can make the land difficult or impossible to see. Cloud cover affects 
the view of the ground and causes problems when using a normalized difference vegetation index (NDVI), 
which in turn hinders the estimation of LST. 


2.2. Ortho-rectification and mosaicking 

Ortho-rectification is a process that removes the geometric distortions introduced during image 
capture and produces an image product that has planimetric geometry in order to precisely register to a ground 
coordinate system and the image scale. After computing geometric models, the final step is to generate 
ortho-rectified strips using a digital elevation model (DEM) and then to mosaic them into a colour-balanced 
image. The information about DEM is also discussed. 


Landsat 8 | 


TIR Data | OLI Data | 
TIR Band 10 & 11 OLI Band 1-9 


Cloud Removal 


Radiance 


Bnghtness Color Balancing Vegetation 
Temperature 
Surface 
Emissivity 
Mono Window Algonthm 
Land Surface Temperature 
Estimation 


Figure 2. Mono window algorithm used to retrieve LST from Landsat 8 imagery 
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2.3. Processing temperature 

Identification of heat anomalies can be performed using analysis of surface temperature according to 
the natural appearance. The Landsat 8 satellite has two sensors, the operational land imager (OLI) sensor and 
thermal infrared sensor (TIRS) sensor with 11 bands. The TIRS sensor provides two thermal bands which 
provide better signal-to-noise (SNR) radiometric performance quantised over a 12-bit dynamic range [17]. 
The OLI sensor contains 9 bands (1-9) and the TIRS sensors have 2 bands (10 and 11). The thermal bands 
(10 and 11) are vital for detecting and analysing the contrasts in land surface temperatures with 100 m 
spatial resolution. 


2.4. Extraction land surface temperature 

The following procedure and formula were used to identify the surface temperature, generate 
the temperature colour map and analyse the data. The thermal band (10.6-11.2 um) (band 10) of the landsat 8 
(OLI) sensor derived the LST of the study area, and the following formula was used to calculate the spectral 
radiance from the digital numbers (DN). The spectral radiance (LA) is calculated with the equation [18]. 


2.5. Top of atmosphere (ToA) spectral radiance 

Digital numbers (DNs) of bands 1-9 were converted to top of atmosphere (ToA) and the spectral 
radiance (LA) was calculated by multiplying multiplicative radiometric rescaling factor of the TIR bands with 
its corresponding TIR band and by adding an additive rescaling factor using (1).These results are then 
compared to DNs reported by the sensor to provide a set of bias and gains (radiometric calibration coefficients) 
for the sensor spectral bands. 
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L=M Q at A: (1) 


where La is the spectral radiance in watts/(m—2 srad—1 um—1); Mz the band-specific multiplicative rescaling 
factor obtained from the metadata (0.000342); A; is the band-specific additive rescaling factor obtained from 
the metadata (0.1); Qca 1s the DN value for the quantized and calibrated standard product pixel of band 10. 


2.6. Brightness temperature (BT) 

BT refers to the amount of electromagnetic radiation that travels upwards from the top of 
the atmosphere [19]. Thermal calibration involves converting the thermal DN values of the raw thermal bands 
from the TIR sensor into TOA spectral radiance, followed by the BT equation shown in (2), 


K2 


TLandsat = mEn) — 273 (2) 
where: Tanaga Effective temperature (Celsius), La Spectral radiance watt/(m**srad*um), K1, K2 Calibration 


constants, € emissivity value given when NDVI process yet. 


2.7. Estimation of emissivity 

Emissivity refers to the radioactive properties of an object and summarises its ability to emit 
radiation [20, 21]. To obtain the LST it is necessary to measure land surface emissivity. For this research, 
the emissivity of the land surface was obtained using NDVI as suggested by Zhang [22]. The corresponding 
land emissivity values as shwon in Table 2 were then calculated from the NDVI results shown in (3). 


Table 2. Estimation of land surface temperature from NDVI [22] 


NDVI Land surface emissivity (€ ) 
NDVI < -0.185 0.995 
-0.185 < NDVI < 0.157 0.970 
0.157 < NDVI < 0.727 1.0094+0.0047 * Ln(NDVI) 
NDVI > 0.727 0.990 


The formula emissivity value given when NDVI process yet: 
_ f: 
E€ = 0.004 `v + 0.986 (3) 


-( NDVI —NDVI,,, ) 
Tom (eh E 


NDVI. — NDVI . 
max min (4) 


2.8. Calculation NDVI 

NDVI is a straightforward graphical indicator that analyses remote sensing measurements and 
assesses whether or not an area contains green vegetation. Studies show that NDVI provides accurate 
estimations of surface radiant temperatures [23-26]. NDVI is one of the outputs that was entered into the model 
to retrieve LST by using (5): 


NDVI = NIK KR 
NIR + R (5) 


where NIR is Near Infra-red band and R is the red band. 


2.9. Land surface temperature (LST) 

LST must be known in order to estimate LST, since it’s a proportionality factor that scales blackbody 
radiance Planck’s law (Ê ) to predict emitted radiance, and it is the efficiency of transmitting thermal energy 
across the surface into the atmosphere [27]. The last step of retrieving the LST or the emissivity corrected land 
surface temperature (celsius) described computed as follows [27]: 
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L Celcius) = T kelvin) 7 272. 15 (6) 
LST = — I (7) 
(1+ws(2)*In(e)) 


where, w = wavelength, ” = 1.438*107 mK 


3. RESULTS AND ANALYSIS 
3.1. Colour balancing 

The Landsat 8 image satellite images have been processed into red (R), green (G), and blue (B) colour 
composites (balancing) with a spatial resolution of 15 m. A color composite of RGB were used as basis of 
visual interpretation. Figure 3 (a) shows the colour combinations. The RGB 432 band (natural colour) denotes 
natural object like vegetation since the colour is green, while RGB 5, 6, 7 showed false color. Both composites 
were used for volcanic products and lineament analyses. And also, to establish the density level of vegetation, 
the red image is useful to gather information about the size of the forest and the maturity of the vegetation as 
showin in Figure 3 (b) covers for estimating the leaf condition related to geothermal system beneath the surface. 





Figure 3. West Sumatra area image processing with colour balancing; 
(a) RGB band 432, (b) RGB band 567 


3.2. SRTM DEM and countour 

Various types of data like shuttle radar topography mission (SRTM) data and digital elevation model 
(DEM) data were used to prepare LST maps of the area. The DEM was derived from topographical maps, and 
the geographical coordinate system and WGS 84 datum were used to perform single map rectification. 
The geo-referenced topographical maps of the area were mosaic. The original SRTM DEM was used for 
the contouring procedure. SRTM data requires a lot of pre-processing to make sure there are no spurious 
artefacts which could become problematic during the analysis, for example, pits, spikes, and patches of no data. 
The final map represents the digital elevation model for all of West Sumatera with 30 m spatial resolution 
as shown in Figure 4. Based on Figure 4 derived from SRTM, some parts of West Sumatera are located directly 
on or extremely close to active volcanoes (Talang, Merapi, Singalong, Sago). Figure 5 shows the contour 
characteristics which were derived from the DEM data. Contours connect points with equal values 
like elevation, temperatures, precipitation, pollution, and atmospheric pressure. The distribution of polylines 
in the contours represent the variation in values over a surface area. 
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Figure 4. SRTM DEM with spatial resolution 30 m 
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Figure 5. Contour characteristics derived from DEM SRTM 


Air temperature decreases when ascending through the troposphere and into the high atmosphere. 
The rate of lapse shows that the temperature falls by 0.650C/100 meters. The fall in temperature at higher 
altitudes is because the earth's surface acts as an energy source that heats the air and makes the temperature 
higher on the planets surface. The land surface temperature estimates that were derived from the satellite 
imagery indicate patches of high and low temperatures, as shown in Figures 6 (a) and 6 (b). 

The calculated land surface temperature as shown in Figures 6 (a) and (b) that in West Sumatera, 
the temperature ranged from -8.6 C° to 32.5 9C°, which is represented by graduated pink to brown shading. 
The areas that were predicted to be higher in temperature are clearly identified (brown colour pixels) in 
the dataset. The Lima Pulu Kota, Tanah Datar, Solok and South Solok city areas show high-temperature values 
(Brown) ranging from 28.1 C° to 32.59 C° because of their high thermal influencing properties. Temperatures 
were lower (Pink) in the north-east region and the range distribution was from -8,6 C° to 5 C°. The hill-shade 
function as shown in Figures 6 (a) and 7 (a) provides the hypothetical illumination of a surface by establishing 
a value for each of the cells contained in a raster. By working out the illumination value of each of the cells 
and setting a hypothetical light source, visualisation of surfaces can be enhanced for analysis or graphical 
display. In Figure 7, show the NDVI also used to identify the vegetation covers for estimating the leaf condition 
related to geothermal system beneath the surface manifestation. 
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Figure 6. (a) Map of LST derived with DEM SRTM 30 m hill-shade, (b) without hill-shade 
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Figure 7. (a) Map of NDVI Overlay DEM SRTM 30 m hill-shade (band 2, 3, 4, 5), (b) without hill-shade 
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4. CONCLUSION 

The thermal imagery from Landsat-8 OLI/TIRS bands can be used to detect surface temperature, 
which in turn allows for the identification of geothermal energy potential depend on surface parameters 
temperature, emissivity and also atmospheric correction. From this method, the brightness temperature 
was derive from LST map of West Sumatra shows that ranged from -8.6 C° to 32.59 C? and the different 
temperatures are represented. The idendtification area from calculated result clearly identifies the hot 
areas in the dataset such as Lima Puluh Kota, Tanah Datar, Solok, and South Solok areas showed 
the high-temperature value in the range of 28.1 C? to 32.59 C° which means that they possess high potential 
for generating thermal energy. The possibility of accurately and efficiently approach to geothermal area 
detection combining with geological survey ananlysis was consideration for future research. 
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